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ABSTRACT 

Detailed in situ data from cirrus clouds have been collected during dedicated field campaigns, but the use of 
the size and habit distribution data has been lagging in the development of more realistic cirrus scattering models. 
In this study, the authors examine the use of in situ cirrus data collected during three field campaigns to develop 
more realistic midlatitude cirrus microphysical models. Data are used from the First International Satellite Cloud 
Climatology Project (ISCCP) Regional Experiment (FIRE)-1 (1986) and FIRE-11 (1991) campaigns and from a 
recent Atmospheric Radiation Measurement (ARM) Program campaign held in March-April of 2000. The 
microphysical models are based on measured vertical distributions of both particle size and particle habit and 
are used to develop new scattering models for a suite of moderate-resolution imaging spectoradiometer (MODIS) 
bands spanning visible, near-infrared, and infrared wavelengths The sensitivity of the resulting scattering prop- 
erties to the underlying assumptions of the assumed particle size and habit distributions are examined. It is 
found that the near-infrared bands are sensitive not only to the discretization of the size distribution but also to 
the assumed habit distribution. In addition, the results indicate that the effective diameter calculated from a 
given size distribution tends to be sensitive to the number of size bins that are used to discretize the data and 
also to the ice-crystal habit distribution. 


1. Introduction 

The general approach for inferring cirrus optical and 
microphysical properties from satellite imagery is to 
compare measured satellite radiances with the results 
of radiative transfer calculations for various conditions 
of viewing geometry, solar illumination, and cloud ma- 
crophysical and microphysical properties. With this ap- 
proach, the retrieved cirrus properties depend on the 
single-scattering properties of the cirrus models used 
in the analyses. In turn, the single-scattering properties 
depend on the assumed cirrus microphysical properties, 
for example, distributions of cirrus particle sizes and 
shapes (or habits). As new cirrus cloud observations 
become available, the cirrus models used in the for- 
ward calculations need periodic reevaluation. The pri- 
mary goal of this study is to document improvements 
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in the development of single-scattering properties for 
a variety of midlatitude cirrus clouds. The cirrus mod- 
els are based upon measurements of midlatitude cirrus 
clouds collected during various field experiments such 
as the First International Satellite Cloud Climatology 
Project Regional Experiment (FIRE): FIRE-I and 
FIRE-II in 1986 and 1991, respectively, and the De- 
partment of Energy (DoE) Atmospheric Radiation 
Measurement (ARM) Program intensive operational 
period (IOP) in 2000. 

The effort to improve the cirrus models is prompted 
in part by recent research discussing the importance of 
cirrus size distribution and particle habits on satellite 
imager-based visible and near-infrared reflectances 
(e.g., Knap et al. 1999; Rolland et al. 2000). New cirrus 
in situ datasets are now available that permit detailed 
examination of the vertical variability of both particle 
size and habit. Another motivation for developing new 
cirrus models is a desire to merge the improved analyses 
of in situ data with improvements in the development 
of single-scattering properties of randomly oriented, 
complex ice-crystal habits such as plates, bullet rosettes, 
aggregates, and columns. 

Recent studies point to the importance of ice-crystal 
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Table 1. Lloud-top and cloud-base temperatures for in situ data 
cases. 


Date 

Cloud-top temperature 
(°C) 

Cloud-base temperature 
CC) 

22 Oct 

FIRE- 1 (1986) 
-37 

-15 

25 Oct 

-34 

-15 

28 Oct 

-37 

-20 

1 Nov 

-44 

-17 


A'i 





25 Nov 

FIRE-II (1991) 
-55 

-32 

26 Nov 

-55 

-22 

5 Dec 

-65 

-39 

9 Mar 

ARM (2000) 
-50 

-30 

13 Mar 

-50 

-20 


habit in retrievals of ice water path, particle size, and 
optical depth in cirrus clouds. The dependence of re- 
trieved optical thickness on crystal habit has been dem- 
onstrated by Mischenko et al. (1996), Doutriaux-Bou- 
cher et al. (2000), and others. One of the most complete 
in situ radiation measurement studies was by Francis et 
al. (1999). They discussed aircraft measurements of 
midlatitude cirrus collected by the Met Office C-130 
aircraft during a flight on 9 November 1995 off the east 
coast of Scotland. In situ microphysical data were ob- 
tained within a thin layer of frontal cirrus, together with 
downwelling filter-band radiometer measurements at 
0.87, 1.61, 3.7, 8.55, and 11 /cm. Banked orbits were 
flown below the cirrus to sample at scattering angles 
ranging from 21° to 127° relative to the sun. To infer 
optical thickness and particle size, radiative transfer cal- 
culations were presented for monodistributions of five 
individual crystal habits: randomized polycrystals 
(Macke et al. 1996), hexagonal columns, hexagonal 
plates, six-branched bullet rosettes, and aggregates com- 
posed of hexagonal columns. For a cirrus cloud com- 
posed of a single habit, single-scattering properties were 
obtained by integrating over size distributions using ef- 
fective radii (defined as the ratio of total volume to total 
projected area) ranging from 9.4 to 85.9 /xm. A ray- 
tracing technique was employed to derive single-scat- 
tering properties of the ice crystals. Francis et al. (1999) 
found that the cirrus composed of polycrystals produced 
effective sizes that were consistent with the in situ data 
across all the wavelengths considered. However, the re- 
trieved 0.87- jim optical thicknesses were considerably 
less than those derived at 1 1 /x m, implying a problem 
with the phase function at solar wavelengths for the 
polycrystal over a wide range of primarily side-scatter- 
ing angles. Francis et al. (1999) explain that the error 
in effective size retrieval introduced by the polycrystal 
overprediction of side scattering is canceled by its un- 
derprediction of single-scattering albedo. Francis et al. 
also note the importance of characterizing realistic cirrus 


Table 2. Spectral and radiometric characteristics of the Earth Ob- 
serving System Terra platform MODIS spectral bands used in this 
study. 


Band No. 

Central 

wavelength 

(#xm) 

Principal absorbing components 

1 

0.65 

H : o, 0 2 0, 

6 

1.64 

H,0, C0 2 CH 4 

7 

2.13 

H,0, CO, CH, N 2 0 

20 

3 78 

HA CO, CH 4 

26 

1.38 

H,0 

29 

8.52 

H,0, O, CH 4 NjO 

31 

11 

H,0, CO, 

32 

12 

H,0, CO, 


ice particles and show that their retrieved optical thick- 
nesses compared better to theory when a laboratory- 
derived phase function reported by Volkovitskiy et al. 
(1980) was employed. 

Baran et al. (1999) used the dual-view geometry of 
alongtrack scanning radiometer (ATSR)-2 data to esti- 
mate the most likely dominating crystal shape at 0.87 
/xm. In their study, radiative transfer (RT) calculations 
were performed for cirrus composed of a single crystal 
type: hexagonal plates, hexagonal columns, six- 
branched bullet rosettes, or randomized polycrystals. 
The simulated radiances subsequently were compared 
to selected ATSR-2 data from a tropical convective cir- 
rus case and a midlatitude cirrus case. For the tropical 
case, the best comparison between the ATSR-2 data at 
0.87 /x m and RT calculations was obtained for the ran- 
domized polycrystal, perhaps due to the lack of en- 
hanced intensity features in the phase function at back- 
scattering angles. For the midlatitude cirrus case, com- 
parison was made between extinction optical depth re- 
trievals using ATSR-2-measured radiances, aircraft 
measurements of midlatitude cirrus at 0.87 and 11 /xm, 
and in situ-calculated optical thickness (Francis et al. 
1999). It is interesting that the best comparison was 
obtained using a phase function obtained from labora- 
tory measurements (Volkovitskiy et al. 1980) rather than 
from the randomized polycrystal. Baran et al. (1999) 
suggested that the aircraft samples the phase function 
at mostly side-scattering angles rather than backscat- 
tering angles, suggesting that perhaps the polycrystal is 
less optimal at side-scattering angles. 

Given the importance of cirrus microphysics in re- 
mote sensing applications, we now turn attention to 
the current set of operational moderate-resolution im- 
aging spectroradiometer (MODIS) cirrus models 
(Baum et al. 2000). The microphysical data are derived 
from the following sources. The cirrostratus and cirrus 
uncinus distributions are described by Heymsfield 
(1975), and Heymsfield and Platt (1984) present two 
modified size distributions typical of warm and cold 
cirrus clouds and three additional cirrus distributions 
that represent midlatitude cirrus at temperatures T = 
— 20°, —40°, and — 60°C. The seven size distributions 
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Table 3 Particle number distributions based on five size bins. 

Max Bin- Number concentration (L~' gnr 1 ) 

dimension width 

(pm) (pm) FIRE-I FIRE-II ARM 

15 20 4 4904 X 10-' 9.2994 X 10~' 0.0000 X 10° 

45 40 1 0089 X 10° 4.2441 x 10" 1 1.0792 X 10° 

130 100 7.6280 X 10“’- 9.8850 x lO" 1 7.2287 x 10" 1 

275 260 2.1286 x 10-’- 1.4585 x 10" : 4.3616 x 10~ 3 

750 640 1.7600 X 10"’ 1.1579 X 10"* 3.5373 X 10~’ 


cover a range of cloud microphysical properties and 
dynamical regimes, for example, production in slowly 
ascending parcels to 100 cm s _1 updrafts. Most of the 
cirrus were sampled aboard a Cessna Citation jet using 
three different probes (Heymsfield 1975; Heymsfield 
and Platt 1984): 1) a forward-scattering spectrometer 
probe (FSSP) for particle sizes ranging between 2 and 
30 /xm in 2-/xm intervals; 2) a one-dimensional optical 
array probe (1D-C) to obtain particle sizes between 20 
and 300 /im in 20- /zm intervals, and 3) a one-dimen- 
sional precipitation probe (1D-P) to size particles be- 
tween 200 and 3000 /im in 200- /am intervals. AForm- 
var (polyvinyl formal — a plastic) replicator was em- 
ployed to provide data on particle habits for crystals 
less than 200 /xm in size; larger particles tended to 
break up upon impact. 

Given the limitations of the 1D-C and 1D-P mea- 
surements, the size spectra for the MODIS cirrus models 
in Baum et al. (2000) are limited to five discrete size 
regions. When the maximum dimension D is small ( D 
< 70 gm), that is, for small ice crystals, the percentage 
of ice crystals within the size region is assumed to be 
composed of 50% bullet rosettes, 25% hexagonal plates, 
and 25% hollow hexagonal columns. When D is greater 
than 70 ^ra, the percentage of ice crystals within the 


size region is assumed to be composed of 30% aggre- 
gates, 30% bullet rosettes, 20% hexagonal plates, and 
20% hollow hexagonal columns. The mixture of habits 
in this initial set of cirrus models is a function only of 
particle size and is invariant between models. 

Recent improvements in the instrumentation used to 
obtain in situ cirrus measurements provide, for the first 
time, a basis for discussing vertical variability of size 
and habit distributions. The cirrus measurements from 
FIRE-II are from a balloonbome ice-crystal replicator 
(Heymsfield and Miloshevich 1995), which provides 
high-quality image, size, and concentration informa- 
tion for particles 10 /xm and above. With this replicator, 
ice crystals have less of a tendency to break up upon 
collection because of the slow rise rates (4-5 m s~') 
of the balloon (Miloshevich and Heymsfield 1997). The 
cirrus measurements from FIRE-I (Heymsfield et al. 
1990; Sassen et al. 1994) and ARM are based on two- 
dimensional crystal cloud (2D-C) and precipitation 
(2D-P) probe data. The 2D-C probe sized particles 
from 25 to about 1000 gm in 25- /xm increments, and 
the 2D-P probe sized particles from 100 to above 2000 
gm in 100-gim increments. The 2D-C and, especially, 
the 2D-P are insufficient to obtain much information 
on habits, especially for the smaller crystals. For FIRE- 
1, habit information was provided through collection 
of particles directly on oil-coated slides (Sassen et al. 
1994). For the ARM campaign, cirrus-particle habit 
data were provided through use of the SPEC, Inc., 
cloud particle imager (CPI), which provides detailed 
imagery of particles, beginning at about 40 gm with 
2.3-gzm resolution. 

With recent advances in developing single-scattering 
properties for ice crystals and the improvements in in 
situ data from cirrus clouds, it is now possible to develop 
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Fig. i. Size distributions based on five size bins for the (a) FIRE-I, (b) FIRE-II, and (c) ARM cirrus microphysical models described in 

section 3. 
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Fig. 2. Size distributions based on 27 size bins for the (a) FIRE-I, (b) F1RE-II, and (c) ARM cirrus microphysical models described in 

section 3. 




more complex cirrus scattering models that include 
more realistic particle size and habit distributions. In 
this study our intent is to form a new set of cirrus models 
based upon the particle size distribution and habit dis- 
tribution observed by in situ measurements of cirrus. 
Previous studies involving in situ cirrus measurements 
have tended to focus on variability in the horizontal, 
whereas we wish to emphasize the vertical variability. 
The new cirrus models represent a significant improve- 
ment in the treatment of both the size distribution and 
the percentages of particle habits. The new cirrus models 
are based on 27 size bins, whereas the initial set of 
models was based on five size bins. The percentage of 
crystal habits for each size bin is based on reanalyses 
of in situ measurements from visual examination of slide 
collections (FIRE-I), balloonborne replicator data 
(FIRE-II), or CPI data (ARM). 

Section 2 presents the data and models used in our 
analyses. Section 3 outlines the development of the new 
cirrus models. Results of RT calculations of top-of-at- 
mosphere reflectances and brightness temperature dif- 
ferences are presented in section 4. Section 5 summa- 
rizes and concludes our study. 

2. Data and models 

a. Balloonborne replicator data: FIRE-II 

distributions 

The FIRE-II size distributions were measured on 25 

and 26 November and 5 December 1991 by the National 
Center for Atmospheric Research (NC AR) balloonborne 
Formvar replicator (Miloshevich and Heymsfield 1997). 
Cloud characteristics for each of the FIRE-II cases are 

listed in Table 1 . The size distribution data are discre- 
tized into 27 size regions ranging from 10 /ira to the 


largest measured size of 540 /am at approximately 100- 
m vertical layers through the cloud, with a visual clas- 
sification of the two dominant crystal habits at each 
level. Images of vertical crystal profiles apparently show 
the presence of three vertical layers (Heymsfield and 
Miloshevich 1995): a production layer composed mostly 
of pristine crystals such as columns or plates, a growth 
layer, and a sublimation layer composed of large and 
often aggregated crystals. 

An average number of particles per size bin is de- 
termined by integration over the height of each cloud 
layer as follows: 

cloud base 

2 n(D, l)Az(l) 

n(D) = ^Z b ~ , (1) 

2 Az(/) 

/= cloud top 

where / is a layer, A z(/) is the thickness of layer /, 
and D is the center maximum dimension of the size 
bin. To determine the layer-averaged percentage of 
habit h at size D, the habit information is averaged 
over height: 

clo ud b ase 

2 f(h, D, l)n(D, l)Az(l) 

m D) = . ( 2 ) 

2 n(D, l)Az(l) 

1 = cloud top 

where f{h, D. I ) is the percentage of habit h at size D 
and layer /. 

b. Aircraft data 

1) FIRE-I DISTRIBUTIONS 

The FIRE-I data were collected from the NCAR King 
Air turboprop in October and November of 1986 over 
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Fig. 3. As in Fig. 2, but the maximum crystal size is limited to 1200 pm so that more detail is presented for the small-crystal size range 


Wisconsin. Ice particle size spectra were derived from 
Particle Measuring Systems, Inc., 2D-C and 2D-P im- 
aging probe data. The 2D-C and 2D-P probes on the 
King Air measured particle size in 25- and 100-gun 
intervals, respectively. A careful examination of slide 
collections provided information on the ice-particle hab- 
it percentages and the presence of small particles. In 
general, few small particles were noted in the collec- 
tions. The data from the FSSP probe (2-30 gm) were 
examined to see whether there were any pockets of liq- 
uid water or whether anomalously high concentrations 
(more than a few tenths per cubic centimeter) were pre- 
sent, which might suggest the presence of relatively high 
concentrations of ice particles. No such high concen- 
trations were detected, suggesting that relatively low 
concentrations of total (and sub-2D-C size) particles 
were present. We have not attempted to account for the 
potential underrepresentation of small particles in the 
2D-C data. 

Cirrus cloud measurements on 22, 25, and 28 October, 
and 1-2 November of 1986 are used in this study. The 
aircraft climbed to cloud top in each case and then con- 
ducted a slow, Lagrangian-type spiral descent down- 
ward to cloud base. All FIRE-I data used were collected 
during the Lagrangian spiral descents. Cloud charac- 
teristics for each of the FIRE-I cases are listed in Table 
1. The merged number concentration datasets use the 
2D-C probe data for particle sizes between 25 and 600 


gm and the 2D-P probe for 600-4000 gun. The number 
concentrations and habit distributions throughout the 
cloud are merged to create an averaged single-layer 
cloud model size distribution by 

clo ud b ase 

2 n(D, j) 

n(D) = (3) 

' V ob, 

and habit distribution 


f(h, D, j) 


f(h. D) = 


respectively, where tV obs is the number of observations 
and j represents a single flight track through the cloud 
layer. 

2) ARM IOP DISTRIBUTIONS 

The ARM data were collected on 9 and 13 March 
2000 by the University of North Dakota Citation aircraft 
near the ARM Cloud and Radiation Test Bed site in 
Oklahoma. Cloud characteristics are listed in Table 1. 
The data from 9 March are of a Lagrangian spiral, and 
those of 13 March are for the entire flight. Ice-particle 
size spectra were derived from 2D-C and 2D-P probe 
data. The merged number concentration dataset uses the 


Table 4. Average maximum dimension and effective diameter based on the microphysica! models described in section 3. 


Field campaign 

FIRE-I 

FIRE-11 

ARM 



5 size-bin fixed habit distribution 


27 size-bin fixed habit distribution 27 size-bin measured habit distribution 


(D) (Atm) 


D, (gm) 


(D) (/tm) 


D, (gun) 


(D) (Atm) 


D, (Atm) 
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Table 5. Particle number distributions for 27 size bins. 


Max 

Bin- 

Number concentration (L 

. 1 pun—) 


dimension widtn 






(pun) 

(pun) 

FIRE-1 

1 FIRE-11 

ARM 


3.0 

4.5 

0.0000 x 

10° 0 0000 x 

10° 

0.0000 x 

10° 

7.5 

6.0 

0.0000 x 

10° 0.0000 x 

10° 

0 0000 x 

10” 

15.0 

8.75 

0.0000 x 

10° 7.1873 x 

10-' 

0 0000 X 

10° 

25.0 

10.0 

8.9809 X 10-' 1.2310 X 

10° 

0.0000 x 

10" 

35.0 

10.0 

1.7962 X 

10° 7.0759 x 

10- 

9.6154 x 

10- 

45.0 

12.5 

1.4813 x 

10° 3.5406 X 

10-' 

7.6923 x 

10- 

60.0 

17.5 

2.2162 x 

10- 3.1284 x 

10- 

1.3679 x 

10“ 

80.0 

200 

1.3541 X 10-’ 1.7818 X 

10- 

2.3871 x 

10” 

100.0 

25.0 

8.6352 X 10- 1.1981 X 

10-' 

6.5310 x 

10- 

130.0 

37.5 

5.5901 X 

10- 2 6.3147 X 

10- 

1.7268 x 

10- 

175.0 

47.5 

3.7979 X 10- 2 5.4743 X 

10- 

9.9477 X 

10- 

225.0 

50.0 

3.0645 X 

10- 2 2.3840 X 

10- 

6.5906 X 

10- 

275.0 

50.0 

^.4290 X 

!0- 2 1.1160 x 

i n—1 
IV 

*♦.4610 a 

nr- 

325.0 

50.0 

1.7150 X 

10- 4.8503 x 

10- 

2.8637 X 

10- 

375.0 

50.0 

1.1570 X 

10- 2.2416 X 

10- 

1.9224 X 

10- 

425.0 

50.0 

7.2273 x 

10-’ 1.5066 x 

10- 

1.4565 X 

10- 

475.0 

62.5 

4.2350 x 

10-’ 5.2346 X 

10- 

1.0156 X 

10- 

550.0 

87.5 

2.4228 X 

10- 1.2870 X 

10- 

7.7074 X 

10- 

650.0 

100.0 

1.7228 x 

10“ 3 0.0000 X 

10° 

3.0205 X 

10- 

750.0 

100.0 

1.7228 x 

10- 0.0000 x 

10° 

1.8922 x 

10- 

850.0 

100.0 

6.9919 X 

10- 0.0000 x 

10° 

1 0486 x 

10- 

9500 

150 0 

5.7010 x 10- 0.0000 X 10° 

4.2422 x 

10- 

1150.0 

225.0 

2-6070 x 

10- 0.0000 x 

10° 

1.8272 X 

10- 

1400.0 

300.0 

1.1388 x 

10- 0.0000 x 

10° 

5.3236 x 

10- 

1750.0 

550.0 

3.4327 x 

10-' 0.0000 X 

10° 

1.0878 x 

10- 

25000 

875.0 

5 3546 x 

10- 0.0000 x 

10° 

1.4242 x 

10- 

3500.0 

1000.0 

2.2661 x 

10- 0 0000 x 

10° 

7.2145 x 

10- 


2D-C probe data for particle sizes between 25 and 900 
pm and the 2D-P probe data for 900-4000 ptm. We 
have not accounted for the potential underrepresentation 
of small crystals in the 2D-C data. Averaged single- 
layer cloud model size distributions are determined from 
Eq. (4). Particle habit distributions were derived from 
visual examination of SPEC, Inc., CPI data. Virtually 
all of the particles were bullet rosettes with an average 
of five bullets per rosette. Larger particles were aggre- 
gates of rosettes. 

c. Single-scattering properties 

An extensive library of single-scattering properties 
has been prepared that encompasses five randomly ori- 
ented crystal habits, 27 size bins, and 10 subbands with- 
in each of the MODIS band numbers in Table 2. The 
crystal habits are solid hexagonal columns, hexagonal 
plates, hollow hexagonal columns, two-dimensional 
bullet rosettes, and moderately rough aggregated col- 
umns. Because the columnar crystals in the King Air 
slide collections and in the replicator data are hollow, 
all columnar crystals are treated as being hollow col- 
umns. The size bins span a maximum dimension range 
of 3-3500 pm. 

The single-scattering calculations are derived using a 
combination of the improved geometric optics method 
(GOM2) and the finite-difference time domain (FDTD) 
techniques discussed in Yang and Liou (1996a.b). The 


Table 6 . Fixed ano measured habit distributions. All particle sizes 
D are in micrometers. 



Fixed 

percentages 

FIRE-1 

FIRE-11 

ARM 

First size limit 

D <10 

fim 

D < 150 

p.m 

D < 100 

D < 100 

fim 

Bullet rosettes 

50% 

37% 

35% 

0% 

Hexagonal plates 

25% 

0% 

46% 

75% 

Hollow columns 

25% 

63% 

17% 

6% 

Aggregates 

0% 

0% 

3% 

19% 

100 pun < 


D > 70 

D > 150 

D > 100 

D < 800 

Second size limit 

/un 

fim 

fim 

fim 

Bullet rosettes 

30% 

33% 

38% 

100% 

Hexaeonal plates 


0% 

o % 

rwc 

u ev 

Hollow columns 

20% 

27% 

22% 

0% 

Aggregates 
Third size limit 

30% 

40% 

40% 

0% 

D > 800 
pun 

Bullet rosettes 

0% 

0% 

0% 

75% 

Hexagonal plates 

0% 

0% 

0% 

0% 

Hollow columns 

0% 

0% 

0% 

0% 

Aggregates 

0% 

0% 

0% 

25% 


FDTD is applied to particles having a size parameter 
smaller than 20, whereas GOM2 is employed for size 
parameters larger than 20. Calculated properties include 
the extinction coefficient, single-scattering albedo, frac- 
tion of delta-transmitted energy, asymmetry parameter, 
and scattering phase function. 

d. Atmospheric absorption 

For this study, we have focused our attention upon 
the creation of correlated ^-distribution routines for the 
spectral regions corresponding to the MODIS channels. 
Because these correlated ^-distribution routines have 
been formulated specifically for the MODIS channels 
rather than being adapted from other studies, we can be 
assured that the routines correctly represent the ob- 
served radiative impact due to molecular absorption 
within these spectral bands. Although our correlated A- 
distribution routines could be applied to an assessment 
of the state of the molecular inventory within the at- 
mosphere, we have used these routines to remove the 
, molecular signal from the data, thereby allowing for an 
analysis of the radiative impacts of other aerosols, 
clouds, surface conditions, and so on. 

A major advantage of using the correlated A-distri- 
bution method is that it can be incorporated directly 
into multiple scattering routines used to consider scat- 
tering as well as absorption by clouds and aerosol par- 
ticles (Lacis and Oinas 1991). Kratz (1995) provides 
an in-depth description of the method used to create 
the correlated A-distribution routines. For both that 
study and our investigation, the derivation of the cor- 
related A distributions has been based upon an expo- 
nential sum fitting of transmissions technique (Wis- 
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wavelengths derived for the cirrus microphysical models described in 
section 3. 




FIRE-1 



FI RE- 11 



ARM 


A (pm) 

5 size-bins 
fixed habit 

27 size-bins 
fixed habit 

27 size-bins 

measured 

habit 

5 size-bins 
fixed habit 

27 size-bins 
fixed habit 

27 size-bins 
measured 
habit 

5 size-bins 
fixed habit 

27 size-bins 
fixed habit 

27 size-bins 
measured 
habit 

0.65 

1 000 

1.000 

1.000 

1.000 

1 000 

1.000 

1,000 

1.000 

1 000 

1.38 

0.9919 

0.9904 

0.9889 

0 9947 

0.9955 

0.9953 

0.9930 

0.9933 

0.9951 

1.64 

0.9063 

0.8834 

0.8835 

0.9334 

0.9418 

0.9406 

0.9182 

0.9214 

0.9399 

2.13 

0.8555 

0.8434 

0.8244 

r\ ont -i 

0071*. 

0.9016 

0.8999 

0.8698 

0.8735 

0.8957 

3.78 

0.6513 

06486 

0.6262 

06834 

0.7062 

0.6995 

0.6632 

0.6777 

0.7113 

8.52 

0.5803 

0.5775 

0.5677 

0.5919 

0.6127 

0.6227 

0.5828 

0.6015 

0.6271 

11.0 

0.5397 

0.5400 

0.5369 

0.5340 

0.5284 

0.5280 

0.5389 

0.5357 

0.5356 

12.0 

0.5680 

0.5682 

0.5648 

0.5653 

0.5635 

0.5641 

0.5677 

0.5671 

0.5694 


combe and Evans 1977) that was applied to the results 
of line-by-line calculations that were run at predeter- 
mined reference pressure and temperature conditions. 
Both studies have relied upon the multiplication trans- 
missivity property to account for the overlap of the 
spectral features associated with the various molecular 
species. We have established the accuracy of the cor- 
related k distribution technique, as compared with the 
line-by-line method, by running atmospheric flux and 
heating rate calculations for a wide variety of atmo- 
spheric conditions. 

e. Radiative transfer model 

Baum et al. (2000) describe the use of a discrete 
ordinates radiative transfer model to calculate the top- 
of-atmosphere (TOA) radiances expected for a range of 
clear and cloudy conditions for the wavelengths listed 
in Table 2. The atmosphere is composed of a discrete 
number of adjacent homogeneous layers. Within each 
layer, the single-scattering albedo and optical thickness 
are constant, but these values may vary from layer to 
layer. For the 1.64-, 1.38-, 2.13-, 3.78-, 8.5-, 11-. and 
12-jU.m channels, the atmosphere is divided into 34 lay- 
ers between the surface and 70 km. For each layer, the 
correlated ^-distribution routines discussed in the pre- 
vious section are implemented to calculate optical thick- 
nesses due to gaseous absorption. For the 0.65-p.m chan- 
nel, for which ozone is the principal absorber and Ray- 
leigh scattering is important, the atmosphere is modeled 
with four layers. The first (uppermost) layer includes 
the effect of ozone absorption and Rayleigh scattering 
from the tropopause to the TOA. The second layer in- 
cludes Rayleigh scattering between the top of the up- 
permost cloud and the tropopause. The third layer con- 
tains the optical properties for either a water or an ice- 
phase cloud. The layer below the cloud contains Ray- 
leigh scatterers. Aerosol contributions are neglected in 
this study. 


for all of the models but the number concentration of 
ice crystals per size bin varies between the models. The 
data upon which we base this study were collected dur- 
ing the FTRE-I (1986), FIRE-II (1991), and ARM IOP 
(2000) field campaigns. As discussed in section 2, the 
data are culled from five cirrus flight missions during 
FIRE-1, three flight missions during hlKE-II, and two 
flight missions during the ARM IOP. 

Two values used to describe a given particle size dis- 
tribution are the average maximum dimension and the 
effective diameter. The average maximum dimension of 
a distribution ( D ) is defined as 

D 

2 Dn(D)A(D) 

(D) = ^ , (5) 

2 n(D)A(D) 

D-D m 


where D mm and D m „ describe the minimum and max- 
imum crystal sizes in the distribution, respectively, and 
A (D) is the width of size bin D. The average maximum 
dimension characterizes the bulk properties of the size 
distribution. The second value is the effective diameter 
D'. Wyser and Yang (1998) demonstrated that the de- 
tails of the size distribution are relatively unimportant 
in specifying the bulk optical properties of cirrus when 
the effective size is calculated as the total volume V 
of the particles divided by the total projected area A. 
This definition is related directly to ice water content 
IWC because V = IW C/p, where p is the density of 
ice. We define D r after Foot (1988) and Francis et al. 
(1994) as 


D 

I 

3 

4 

2 V(h, D)f {h, D)n(D)MD) 

h=l 

2 

2 

4 

2 A(h, D)f (h, D)n(D)MD) 
/»= 1 


3. Development of cirrus models 

a. Cirrus microphysical models 

In the MODIS version- 1 cirrus models described by 
Baum et al. (2000), the habit distribution is the same 


where /(/i. D) is the percentage of a specific particle 
habit for size D, and n(D ) is the number concentration 
of crystals. 

To investigate the effects of habit distribution and 
number of size bins on the cirrus single-scattering prop- 
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Fig. 4 Scattering phase functions based on F1RE-I data for the three cirrus microphysical models described in section 3 for four 

wavelengths: 0.65. 1 64. 3.78. and 11.0 /xm. 


erties (SSP), we created three sets of cirrus models using 
the SSP library. The three model sets are as follows. 

1) The cirrus models, based on FIRE and ARM in situ 
measurements, have five size bins and the MOD1S 
version- 1 fixed habit percentages (i.e., the same habit 
percentages are used for each model). The maximum 
dimension does not exceed 1100 /am. 

2) The cirrus models have expanded size distributions 
(27 size bins) based on in situ data but have the 
MODIS version- 1 fixed habit distribution. 

3) The cirrus models have expanded size distributions 
containing 27 size bins and habit distributions based 
on in situ data. 

These models will now be discussed in detail. 

1) Cirrus models: Five size bins and fixed 

HABIT DISTRIBUTION 

As. a first step in developing cirrus models from in 
situ data, the first set of models is patterned after those 
presented in Baum et al. (2000). The flight-track-av- 


eraged size distributions are discretized into five size 
bins as shown in Fig. 1 and tabulated in Table 3. As 
with the Baum et al. (2000) models, however, the habit 
distributions for each bin are kept constant for each 
model. When the maximum dimension D is small (D 
< 70 /im), that is, for the smallest two size bins, the 
percentage of ice crystals within the size region is 
assumed to be 50% bullet rosettes, 25% hexagonal 
plates, and 25% hollow hexagonal columns. When D 
is greater than 70 /im (i.e., for the largest three size 
bins), the percentage of ice crystals within the size 
region is 30% aggregates, 30% bullet rosettes, 20% 
hexagonal plates, and 20% hollow hexagonal col- 
umns. The D r and ( D ) values for these distributions 
are listed in Table 4. 

2) Cirrus models: Twenty-seven size bins and 

FIXED HABIT PERCENTAGES 

In this section, the in situ data used in the previous 
section are discretized into 27 size bins while the habit 
distribution is left unchanged. The 27 size bins cover a 
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Scattering Angle (°) Scattering Angle (°) 

Fig. 5. Same as Fig. 4 but based on FIRE-II data. 


size range from 3 to 3500 /am. The average number 
concentration n(D) for each size bin is determined by 

.V 

2 n(D, case)A(D) 

n{D) = , (7) 

2 MD) 

casc= 1 

where N refers to the number of individual cirrus cases. 
For FIRE-I and .ARM. the in situ measurements were 
obtained by aircraft, but balloon replicators were used 
in FIRE-II. Where aircraft were used, a case would rep- 
resent the data collected for all the flight tracks for a 
cirrus mission on any given day. When replicators were 
used, a case would be an entire dataset collected for the 
path taken through the cirrus. 

Figure 2 shows the size distributions calculated 
from averages of five cases from FIRE-I, three cases 
from FIRE-II, and two cases from ARM. Figure 3 
shows the same size distributions as in Fig. 2, but the 
size range is limited to 1200 /am to bring out the detail 
of the small size bins. Note that there are no large 
particles recorded by the replicator data as there are 


from the aircraft probe measurements. Table 5 pro- 
vides the number concentrations for the size distri- 
butions. Based on an average of the five cases from 
FIRE-I, the sizes range from 25 to 3500 /am. The 
particle sizes for the FIRE-II replicator data range 
from 15 to 600 /am, and the sizes from the two ARM 
IOP cases range from 35 to 3500 /am. Table 4 provides 
average values of (D) and D t for these cases. 

3) Cirrus models: Twenty-seven size bins and 

MEASURED HABIT PERCENTAGES 

At this point, the in situ data for the cirrus models 
have been discretized into 27 size bins. In this section, 
the in situ data are used to provide averages of habits 
based on particle size as follows. The average distri- 
bution habit weights /(£>) are determined by 

,V 

2 f(D, cas e)rt(D, case)A (D) 

f(D) = ^—7 • (8) 

2 n(D, case)A(D) 

case= 1 

w here N refers to the total number of cases. 
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Fig. 6. Same as Fig 4 but based on ARM data. 


The averaged habit distributions based on the in situ 
measurements is provided in Table 6. We note that al- 
though aggregated crystals are found in situ only for 
crystal sizes greater than about 300 gm, we include 
them in our models at smaller sizes to account for par- 
ticles that have complex or rough shapes, even though 
the observed crystals may not technically have been 
aggregates. Table 4 shows the (D) and D r values cal- 
culated when both particle size and habit are incorpo- 
rated into the calculations. 

4) Comparison of model effective diameters 

Comparison of D r values in Table 4 shows clearly 
the sensitivity of both size distribution and habit in the 
calculations. For the FIRE-I data, an increase from five 
size bins (103.1 jam) to 27 size bins (124.2 gm) results 
in an increase of D € by 20%, whereas a further modi- 
fication in the habit distribution results in an additional 
increase in D, of 1 1%. Much of this difference probably 
can be attributed to the greater range of sizes in the 27- 


bin case (3-4000 gm) than in the five-bin case (5-1070 
gim). The largest crystals are not included in the five- 
size bin distribution. 

For the FIRE-II case, D e decreases by 18% when 
the number of size bins increases from 5 to 27, keep- 
ing the habit distribution fixed. Furthermore, modi- 
fying the distribution of habits has little appreciable 
effect on D t . For the ARM data, an increase the num- 
ber of size bins has almost no effect on D r , changing 
the value by only 1%. However, for the ARM cirrus 
model, the D c shows a surprising sensitivity to particle 
shape. It is important to note that the ARM data are 
composed of measurements from cirrus uncinus that 
were composed primarily of bullet rosettes. Although 
the maximum dimension of a bullet rosette may be 
very large, the ratio of total volume to projected area 
is typically much smaller, resulting in a smaller value 
for D e . 

These results suggest that the retrieval of D, from 
satellite data will depend critically on the assumed dis- 
tributions of both particle size and habit and point to 
the need for more in situ data. 
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b. Cirrus single-scattering properties 

A MODIS band-averaged size distribution is cre- 
ated by integrating single-scattering properties over 
the spectral response function (SRF), habit distribu- 
tion, and size distribution. The SRF weights / A are 


normalized with respect to the discrete spectral sub- 
bands sb such that 

10 

X A(sb) =1. (9) 

sb= 1 

The mean extinction and scattering cross sections are 
given by 


fiw I 

X 

r « 

X fih. D)n(D)MD) 

X aJLK D, sb)/ A (sb) ■ 

sb 

J 

4 

X m D)n(D)A(D) 



rx _ r> Il-i I 

L" _ ' J 


Dun 

X 

D=Di 

r ^ 

X fQi. D)n(D)A(D) 

l*" 

X D, sb)/ A (sb) | 

Dun 

X 

d=d m 

4 

X f(h, D)n(D)A(D) 



( 10 ) 


(ID 


respectively, where h refers to the habit (summed over 
hexagonal plates, hollow columns, bullet rosettes, and ag- 
gregates), / A refers to the SRF weights, and A (D) is the 
width of size bin D. The extinction coefficient is given by 


/3„ 


X 




X f(h, D)n(D)MD) 

h= 1 


The dimensionless optical thickness of the cloud may 
be obtained from this relationship: 

T = P'^Az, (13) 

where A z is the cloud geometrical thickness. 

The averaged single-scattering albedo, delta-trans- 
mitted energy, and phase function at scattering angle 0 
are given, respectively, by 


X 


X D, sb)A(sb) . 


(12) 


cr , 


(14) 


/» 


D^ 

X 

D= D mm 

[x /(*. D)niD)A(D) 

L * 

X f»{h, D, sb)t7- sc .(ft, D, sb)/ A (sb) | 

D 

D=D m ± 

[X /(*. D)n(D)A(D) 

L * 

X <r m (h, D, sb)/ A (sb) 

sb 

1 


and 


P(Q) = 


Cm., 

X 

n=D„- I h 


X x /(*». D)n(D)A(D) 


T P(Q, h, D, sb )trJLK D, sb)/ A (sb) 




X X fOl D)n(D)A(D) 


D=D, n- hub 


X D, sb)/ A (sb) 


(15) 


(16) 
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Table g. Extinction coefficients 0, (km 1 ) for a set of MODIS wavelengths derived for the cirrus raicrophysica! models described in 

section 3 




F1RE-I 



FIRE-II 



ARM 


A (/am) 

5 size-bins 
fixed habit 

27 size-bins 
fixed habit 

27 size-bins 
measured 
habit 

5 size-bins 
fixed habit 

27 size-bins 
fixed habit 

27 size-bins 
measured 
habit 

5 size-bins 
fixed habit 

27 size-bins 
fixed habit 

27 size-bins 
measured 
habit 

065 

0.752 

0.735 

IBI3 ; ' 

rmsm 

0.276 

0.269 

2.169 

1.546 

1.326 

1,38 

0.759 

0.738 

gjgSj- 1 ( 


0.276 


2.168 

1.551 

1.333 

1 64 

0.750 

0.736 


0.362 

0.276 

0.269 

2.175 

1.548 

1.321 

2.13 

0.765 

0.739 


0.369 

0.276 

0.270 

2.196 

1.545 

1.303 

3.78 

0.745 

0.733 


0.358 

0.277 

0.268 

2.153 

1.559 

1.355 

8.52 

0.757 

0.738 


0.357 

0.277 

0.278 

2.169 

1.590 

1.404 

11.0 

0.716 

0.701 


0.338 

0.253 

0.246 

2.065 

1.456 

1.250 

12.0 

0.739 

0.723 

liB 

0.353 

0.268 

0.261 

2.132 

1.516 

1.310 


Table 7 provides the single-scattering albedos w c that 
have been computed for a set of MODIS wavelengths 
for the cirrus microphysical models developed in section 
3. For all cirrus microphysical models, the value of 
is 1 at 0.65 jam, but that value tends to decrease with 
increasing wavelength. At the near-infrared (NIR) wave- 
length of 1.38 jam, the calculations show that w„ does 




Fig. 7. Scattering phase functions for FIRE-1. FIRE-11, and ARM 
size bins, measured habit distributions) for ft 


not depend very strongly on either the discretization of 
the size distribution or variations in the crystal habit 
distributions. 

At NIR wavelengths of 1.64, 2.13, and 3.78 /am, 
where the imaginary index of refraction for ice has a 
much higher value than at 0.65 or 1 .38 /am, the co a values 
are shown to have more dependence on the details of 




a for the cirrus microphysical models described in section 3c (27 
wavelengths- 0.65, 1.64. 3.78, and 1 1.0 /am. 






March 2002 


N A S 1 R I E T A L . 


209 



220 230 240 250 260 270 280 290 300 220 230 240 250 260 270 280 290 300 

1 1 urn BT (K) _e_ 27 size bins, measured habit 11 tim BT (K) 

A 27 size bins, fixed habit 
-a- 5 size bins, fixed habit 

Fig 8, Radiative transfer results generated using the FIRE-I cirrus microphysical models derived in section 3 for (a) 0.65-/im reflectance 
as a function of 11.0- /cm brightness temperature, (b) 1 .64- ^.m reflectance as a function of 11.0-/xm brightness temperature and 0.65-/rm 
reflectance as a function of 1 1.0-^m brightness temperature, (c) BTD (3.78-1 1) as a function of 1 1 0-fim brightness temperature, and (d) 
BTD (8 ,52—1 1) as a function of 1 1 0-^im brightness temperature The RT calculations are performed assuming a cirrus cloud at T = 230 K 
using a climatological midlatitude summer temperature and humidity profile to calculate atmospheric absorption, given solar zenith, viewing 
zenith, and relative azimuth angles of 30°, 20°, and 140°, respectively. 

both the size and habit distributions than at the shorter /am, the w„ values decrease when the habit distribution 

wavelengths. At 1.64 |tm, all three models show an is changed from the fixed distribution to the measured 

increase in u> 0 when the number of size bins is increased distributions for the FIRE-I and FIRE-II models, where- 
from 5 to 27. This is also the trend for the FIRE-II and as the values for the ARM model increase. The FIRE- 
ARM cases at 2.13 and 3.78 /rm. although the trend is II models tend to be slightly more sensitive to the dis- 

reversed for the FIRE-1 case. At 1.64, 2.13, and 3.78 cretization of the size bins, and the ARM and FIRE-I 
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Fig. 9. Scattering angles computed for the RT calculations 
provided in Figs. 10, 12, and 14. 


models are more sensitive to the details of the habit 
distribution. These sensitivities can be explained to a 
large degree by the dependence of w c on D,. 

At IR wavelengths, where ice absorption within the 
crystals is much stronger than at visible or NIR wave- 
lengths, the values for still show sensitivity both to 
the assumed ice crystal size and to the habit distribu- 
tions. The behavior of the a> 0 values with respect to 
changes in wavelength, habit distribution, and number 
of size bins is similar for the FIRE-II and ARM models 
but differs for the FIRE-I model. The sensitivity to the 
discretization of the size bins and to the habit distri- 
bution at IR wavelengths is similar to that at the NIR 
wavelengths. The FIRE-I and ARM modeled io„ values 
are more sensitive to the size distributions; the values 
for the FIRE-II model are more sensitive to the habit 
distribution. 

Scattering phase functions for the FIRE-I, FIRE-II, 
and ARM cirrus microphysical models are provided in 
Figs. 4 through 6, respectively, for MODIS wavelengths 
A = 0.65, 1.64, 3.78, and 11.0 /xm. For the FIRE-I and 
FIRE-II models, there tends to be little variation at scat- 
tering angles greater than 50° among the three sets of 
microphysical models, with the largest variations oc- 
curring at 3.78 jzm. The phase functions developed for 
the ARM data tend to display the most variation, pri- 
marily due to the prominence of bullet rosettes as the 
dominant habit and also the prevalence of very large 
ice crystals. 

A comparison of phase functions developed from the 
use of measured ice crystal habit distributions and a 27- 
bin size distribution is provided in Fig. 7. This com- 
parison shows that the most variation in the phase func- 
tion occurs at 3.78 jr,m, especially at side-scattering an- 
gles. 


4. Radiative transfer results 

Radiative transfer calculations were performed as in 
Baum et al. (2000) to demonstrate the effects of habit 
and size distribution on TOA reflectances and brightness 
temperatures. Results from using a midlatitude summer 
temperature and humidity profile to compute atmo- 
spheric absorption as described in section 2d are shown. 
Although not shown, similar results are obtained when 
a midiatitude winter profile is used. 

Reflectance and brightness temperature difference re- 
sults, based on FIRE-I cirrus data, are provided in Fig. 
8 for solar zenith, viewing zenith, and relative azimuth 
angles of 30°, 20°, and 140°, respectively. The radiative 

trancfpr ralrinltitinnc orp nprfnrmpH fnr o ptmic of 

UMUUlVt VMVW1UWWUO UX W ^/Vi iVi UiVU XV/1 U bill UJ VIOUU Ul 

T = 230 K. Cirrus 0.65-^.m optical thicknesses vary 
from 0 to 20. Figure 8a shows that the 0.65-ptm reflec- 
tances are not a strong function of either size or habit 
distributions. However, the 1. 64- ^im reflectances shown 
in Fig. 8b indicate a dependence both on size and habit 
distribution, and the reflectance differences between the 
models increase with increasing optical thickness. Figs. 
8c and 8d show brightness temperature difference 
(BTD) results between the 3.78- and 1 1-^m bands [BTD 
(3.78-11)] and the 8.5- and ll-/zm bands [BTD (8.5- 
11)]. Of interest, the BTD (3.78-11) calculations show 
more sensitivity to the particle habit and size distribu- 
tions than do the BTD (8.5-11) results. 

The results in Fig. 8 are performed for a specific set 
of viewing angles and will now be extended to a wider 
range of viewing angles. As a reference for the follow- 
ing figures, a set of scattering angles is provided in Fig. 
9. Figure 10 shows the result of RT calculations per- 
formed assuming a cirrus 0.65-/im optical thickness of 
1 at a constant solar zenith angle of 30° for the three 
cirrus microphysical models (see section 3a) derived 
from the FIRE-I data. The numerical differences of re- 
flectances or brightness temperature differences become 
more pronounced at higher cirrus optical thicknesses but 
are not shown here. 

Figure 10a shows the percent relative differences in 
1.64-/xm reflectances between the FIRE-I cirrus model 
derived using 27 size bins with a measured distribution 
[section 3a(3)] and the cirrus model derived using 27 
size bins with fixed habit distribution [section 3a(2)]. 
Figure 10b shows the percent relative differences in 
1.64-/xm reflectances for the FIRE-I cirrus models de- 
rived using 27 size bins with a fixed distribution [section 
3a(2)] and the cirrus model derived using five size bins 
with fixed habit distribution [section 3a( 1 )]. The dif- 
ferences in 1.64-/um reflectance shown in Fig. 10a are 
due primarily to a change in the prescribed habit dis- 
tribution, whereas the differences in Fig. 1 0b are due 
to changing the discretization of the size distribution. 
The reflectance differences in Fig. 10a and 10b can be- 
come appreciable at scattering angles greater than 160° 
and at scattering angles less than 120°, even though the 
cirrus optical thickness is constant at 1. As the cirrus 
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e) 8.5 urn -11 jun f) 8.5 urn -11 pm 



0 30 60 90 120 150 180 0 30 60 90 120 150 180 

Relative Azimuth Angle (°) Relative Azimuth Angle ( c ) 

Fig. 10. Relative difference computed using the FIRE-I cirrus microphysical models. Calculations assume a 0.65 -/cm optical thickness of 
1, a solar zenith angle of 30°, and a midlatitude summer atmospheric profile. 


optical thickness increases. Fig. 8 indicates that the 
1.64- /cm reflectances will show even more sensitivity 
to the prescribed particle and habit distributions. 

Figures 10c and lOe show the differences (in kelvin) 
in the BTD (3.78-11) and BTD (8.5-11) values cal- 
culated for the FIRE-I cirrus model derived using 27 
size bins with a measured habit distribution [section 


3a(3)] and the cirrus model derived using 27 size bins 
with fixed habit distribution [section 3a(2)]. The dif- 
ferences in BTD (3.78-11) and BTD (8.5-11) are rel- 
atively small over the prescribed range in scattering an- 
gles. Examination of Fig. 8 supports the contention that 
the BTD (3.78-1 1 ) differences are expected to become 
more pronounced at higher cirrus optical thicknesses. 
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220 230 240 250 260 270 280 290 300 220 230 240 250 260 270 280 290 300 



220 230 240 250 260 270 280 290 300 220 230 240 250 260 270 280 290 300 

11 pm BT (K) . „ . .. . . ... 11pmBT(K) 

-©- 27 size bins, measured habit p ' 

A 27 size bins, fixed habit 

-o- 5 size bins, fixed habit 

Fig. 11. As in Fig. 8, but using cirrus microphysical models based on FIRE-11 data. 


Figures lOd and lOf show the differences (in kelvin) typically less than 0.8 K at viewing zenith angles less 
in the BTD (3.78-11) and BTD (8.5-11) values cal- than 40°. The BTD (8.5-11) differences are typically 
culated for the FIRE-I cirrus model derived using 27 less than 0.1 K at viewing zenith angles less than 40°. 

size bins with a fixed habit distribution [section 3a(2)] Figure 1 1 shows results similar to those in Fig. 9 but 
and the cirrus model derived using five size bins with for the FIRE-II data. The calculations are performed for 
fixed habit distribution [section 3a(l )]. Here the differ- a cirrus cloud at T = 230 K using a climatological 
ences in BTD (3.78-1 1) and BTD (8.5-1 1) are due pri- midlatitude summer profile. As with the FIRE-I results, 
marily to a change in the discretization of the particle Fig. 11a shows that the 0.65-/xm reflectances are not a 
size distribution. The BTD (3.78-11) differences are strong function of either the assumed size or habit dis- 
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e) 8.5 pm - 1 1 urn f) 8.5 pm -11 pin 



0 30 60 90 120 150 180 0 30 60 90 120 150 180 

Relative Azimuth Angle (') Relative Azimuth Angle f) 

Fig. 12. Relative difference computed using the FIRE-II cirrus microphysical models. Calculations assume a 0 65- rim optical thickness of 
I. a solar zenith angle of 30°, and a midlatitude summer atmospheric profile 


1 l-/um brightness temperature. The BTD (3.78-1 1) cal- 
culations show more sensitivity to the particle size dis- 
tribution than to the habit distribution. The BTD (8.5- 
1 1) results indicate less sensitivity to changes in size or 
habit distribution than with BTD (3.78-11). 

Figure 12 shows similar results to Fig. 10 but for 
the FIRE-II data. Figures 12a, 12c, and 12e show the 


tribution. However, the 1.64-p.m reflectances shown in 
Fig. lib indicate a slight dependence both on size and 
habit distribution at higher cloud optical thicknesses. 
The effect is less pronounced for the FXRE-II data than 
for the FIRE-I data. Figures 11c and lid show bright- 
ness temperature difference results for BTD (3.78-11) 
and the BTD (8.5-1 1). respectivelv, as a function of the 
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A 27 size bins, fixed habit 
-o- 5 size bins, fixed habit 


Fig. 13. As in Fig 8, but using cirrus microphysical models based on ARM IOP data. 


differences between the two cirrus models that each 
have 27 size bins but different habit distributions, 
whereas Figs. 12b, 1 2d, and 12f show the differences 
between the cirrus models that have the same habit 
distribution but a different number of size bins de- 
scribing their size distributions. The 1.64-/xm reflec- 
tances are not strongly influenced by changes in the 
habit distribution over a range of scattering angles, 
with reflectances generally within 3% of each other. 
The BTD (3.78-11) show slightly more dependence 


on the size distribution than on the habit distribution. 
For the FIRE-II data, the BTD (3.78-11) values are 
more sensitive to both size and habit distribution than 
are the BTD(8.52-11) values. 

Figure 13 shows results similar to those in Figs. 8 and 
1 1 but for the ARM IOP data. The calculations are per- 
formed for a cirrus cloud at T — 230 K using a clima- 
tological midlatitude summer profile. As with the FIRE- 
I and FIRE-II results. Fig. 13a shows that the 0.65-yu.m 
reflectances are not a strong function of either the as- 
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Fig. 14. Relative difference computed using the ARM cirrus microphysical models. Calculations assume a 0 65-/xm optical thickness of 
1. a solar zenith angle of 30°. and a midlatitude summer atmospheric profile 


sumed size or habit distribution. The 1.64-jzm reflec- 
tances shown in Fig. 13b indicate a slight dependence 
both on size and habit distribution at higher cloud optical 
thicknesses. Figures 13c and 13d show results for BTD 
(3.78-11) and BTD (8.5-11), respectively, as a function 
of the ll-/xm brightness temperature. The BTD (3.78- 
1 1 ) calculations show more sensitivity to the particle hab- 


it distribution than to the size distribution. The main fea- 
ture of the ARM data is that the predominant habit was 
bullet rosettes at all levels in the cloud. The bullet rosettes 
have a relatively low volume-to-area ratio when com- 
pared with other crystal habits such as plates or columns. 
The BTD (8.5-11) results indicate little sensitivity to 
changes in size or habit distribution. The most sensitivity 
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is shown iii a relatively uaiiow i ante uf optical thick- 
nesses ranging from 0.5 to 3. 

Figure 14 shows results similar to those in Figs. 10 
and 12 but for the ARM IOP data. The 1.64-p.m re- 
flectances show markedly more sensitivity to the change 
in habit distribution (Fig. 14a) than to the size distri- 
bution (Fig. 14b). For the ARM IOP data, the 1.64- /xm 
reflectances are influenced more by the size distribution 
than by changes in the habit distribution. The effect of 
using measured habit distributions (Figs. 14c,e) rather 
than the fixed habit distribution tends to lower both the 
BTD (3.78-11) and the BTD (8.5-11) values. The op- 
posite tendency is noted in Figs. 14d and 14f when the 
number of size bins increases from 5 to 27 while keeping 


5. Summary and conclusions 

Detailed in situ data from cirrus clouds have been 
collected during dedicated field campaigns, but the use 
of the size and habit distribution data has been lagging 
in the development of more realistic cirrus scattering 
models. In this study, we examine the use of in situ 
cirrus data collected during three field campaigns to 
develop more complex and, it is hoped, more realistic 
cirrus microphysical models. Data are examined from 
two FIRE campaigns (held in 1986 and 1991) and a 
recent ARM IOP campaign (held in 1999). All three 
field campaigns were held in the central United States. 

The new cirrus microphysical models form the basis 
for the development of bulk scattering properties such 
as single-scatter albedo, extinction coefficients, and 
scattering phase functions for wavelengths correspond- 
ing to a set of MODIS visible, near-infrared, and infrared 
bands ranging from 0.65 to 12 /am. The sensitivity of 
the resulting scattering properties to the underlying as- 
sumptions of the assumed particle size and habit dis- 
tributions are examined. Our results indicate that the 
effective diameter calculated from a given size distri- 
bution tends to be sensitive to the number of size bins 
that are used to discretize the data and also that the ice- 
crystal habit distribution can play an important role. 

To derive a set of midlatitude cirrus models, our ap- 
proach is to reexamine in situ data from a variety of 
midlatitude cirrus cases for vertical distributions of par- 
ticle size and particle habit. Although this is a time- 
consuming effort, it offers a firm basis upon which to 
derive scattering properties. Several questions arise as 
to our approach. How many individual cirrus cases are 
necessary to obtain a representati ve sample for each type 
of cirrus? Does it make more sense to use cloud-top 
temperature as a defining factor to define a set of cirrus 
models or the origin of the citrus? If one were to adopt 
a certain cirrus microphysical model upon which to de- 
rive a set of scattering properties over a range of wave- 
lengths, would the retrieved cirrus properties be con- 
sistent no matter which wavelengths were used to infer 
the properties? Because we did not try to account for 


the potential underrepresentation of small crystals by 
the 2D-C probes for the FIRE-I and ARM cases, the 
extinction coefficient may be underestimated. This pos- 
sibility raises questions regarding the effects of includ- 
ing parameterizations of small crystals in the size dis- 
tributions. 

One way to test the consistency of results using dif- 
ferent satellite imagers would be to derive cirrus scat- 
tering models in a similar fashion as that presented in 
this study for other satellite imagers such as the Geo- 
stationary Operational Environmental Satellite (GOES), 
Advanced Very High Resolution Radiometer (AVHRR), 
the Multiangle Imaging Spectroradiometer, and the Geo- 
stationary Imaging Fourier Transfer Spectrometer, 
which will provide infrared wavelength remote sensing 
of cirrus. For pairs of these satellites, such as GOES 
and AVHRR or GOES and MODIS, one can retrieve 
cirrus properties for the scene where the data for each 
imager are taken within minutes of each other. If the 
retrieved cirrus properties are consistent across the 
wavelength spectrum offered by these and other instru- 
ments, it would provide some evidence that a set of 
cirrus models derived in such a fashion is representative 
of midlatitude cirrus. 

As a result of the research presented in this study, we 
are now reanalyzing cirrus in situ data from other field 
campaigns to build microphysical and radiative models 
for convective systems, tropical cirrus, and polar cirrus. 
Further efforts will be directed toward the development 
of a suite of cirrus models that can be applied appro- 
priately to any type of ice-cloud system for global sat- 
ellite data analysis. 
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